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Abstract. We present the XMM-Newton first light image, taken in January 2000 with the EPIC pn camera during 
the instrument's commissioning phase, when XMM-Newton was pointing towards the Large Magellanic Cloud 
(LMC). The field is rich in different kinds of X-ray sources: point sources, supernova remnants (SNRs) and diffuse 
X-ray emission from LMC interstellar gas. The observations are of unprecedented sensitivity, reaching a few 10 32 
erg/s for point sources in the LMC. We describe how these data sets were analysed and discuss some of the 
spectroscopic results. For the SNR N157B the power law spectrum is clearly steeper than previously determined 
from ROSAT and ASCA data. The existence of a significant thermal component is evident and suggests that 
N157B is not a Crab-like but a composite SNR. Most puzzling is the spectrum of the LMC hot interstellar 
medium, which indicates a significant overabundance of Ne and Mg of a few times solar. 
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1. Introduction 

On December 10 1999 the XMM-Newton spacecraft was 
placed in a 48 hour Earth orbit by the first commercial 
ARIANE V launcher. XMM-Newton, or X-ray Multi- 



Mirror Mission (Jansen et al. 2001 ), is the second corner- 
stone of the Horizon 2000 science program of the European 
Space Agency ESA. 

The EPIC pn camera (Striider et al. |2001| ) was suc- 
cessfully commissioned during the period of mid January 
through mid March (Briel et al. |2000[ ). On January 19 
2000, shortly after the switch-on of the EPIC pn camera, 
XMM-Newton received first light on the EPIC pn CCD 
chip. We selected as first light target a region in the LMC 
about 10 arcmin southwest of the 30 Doradus A complex. 
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The field has been known to be very rich in X-ray sources 
from Einstein ( Long e t al. 1981 ) and ROSAT observations 
(Triimper et al. 1991 ). There are at least three SNRs (one 
of those known to contain a pulsar) , one OB association 
and two R associations as well as the SN 1987A. This 
XMM-Newton first light has already stimulated detailed 
studies of some of the objects. Here we report some results 
of the observation with four observation sequences of the 
same field added totaling an exposure time of 106 ks. The 
observations are of unprecedented sensitivity, reaching a 
few 10 32 erg/s for point sources. We present and discuss 
the spectroscopic results over the broadest X-ray band- 
pass of 0.3- 12 keV ever covered by a single imaging X-ray 
instrument for the entire field of 27 x 27 arcmin 2 . In section 
2 we describe details of the data analysis techniques which 
may be useful for other observations as well. In section 3 
we discuss some of the scientific results. 
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Table 1. Journal of Observations 



Date 


Time 


RA(2000) 


DEC(2000) 


2000 Jan f9 


16:19 


- 17:23 


05 


36 


57.0 


-69 


13 


47 


2000 Jan 19/20 


17:30 


- 04:28 


05 


36 


57.0 


-69 


13 


47 


2000 Jan 21 


15:38 


- 19:38 


05 


37 


04.0 


-69 


13 


00 


2000 Jan 21/22 


20:32 


- 07:01 


05 


37 


04.0 


-69 


13 


Of) 


2000 Jan 22 


09:08 


- 12:02 


05 


36 


04.0 


-69 


13 


00 



2. Observations and Data Reduction 

The region in the LMC was observed with the pn CCD 
camera operating in full frame mode (frame time: 73 ms), 
with the medium and thin filter inserted during the first 
thr ee an d last two observa tions, respectively (Striider et 
al. |2001| and Turner et al. |200l| ). In Tab. we show the 
journal of the observations. For a fixed pointing direc- 
tion some photons fall onto the gaps between adjacent 
CCD chips. Fortunately, the pointing direction and the 
roll angle were changed for the later three observations 
with about the same exposure time as for the previous 
two observations, such that the surface brightnesses are 
within a factor of 2, with any noise contribution negligi- 
ble. The data analysis was performed with software which 
was developed for the analysis of ground calibration mea- 
surements of the EPIC pn camera (Dennerl et al. 1999| ). 
Most of this software has been implemented into t he XM M 
Standard Analysis Software (SAS, Watson et al. 2001 ). 

During the observations the background induced by 
low energy protons changed and reached sometimes count 
rates of several hundred counts per second in the full 
energy band. To enhance the signal to noise ratio, we 
screened the observations for times with high background 
and excluded those times from further processing. The net 
observing time after screening is 54.8 ks. In addition the 
trails of minimum ionizing particles and reemission trails 
were removed frame by frame. The remaining events were 
corrected for charge transfer loss and gain variations, and 
events created by a single photon, but split across several 
pixels, were recombined, to yield a photon event file in 
detec tor co ordinates with corrected energies (see Dennerl 
et al. 1999 for a detailed description). 

Test images produced from this data set in various en- 
ergy bands indicated that the energy range E = 0.3 — 
5.0 keV was optimal to provide clean images with high 
signal to noise ratio. The energy range was then split into 
three narrower bands of similar signal to noise ratio, to cre- 
ate red, green and blue components of the final rgb image. 
The variation in spectral hardness of the LMC sources and 
the spectral resolution of the EPIC pn camera resulted in 
colourful images which, however, do not represent the dy- 
namic range in intensity very well. To preserve the flux 
information as well, the colour saturation was reduced by 
choosing overlapping energy ranges (0.3-1.0 keV for red, 
0.8-2.0 keV for green, and 1.5-5.0 keV for blue), and by 
further distributing some of the flux from these bands into 
other bands. 




Fig. 1. Superposition of the photons detected with the 
EPIC pn camera in all five LMC observations, displayed 
in equatorial coordinates. The photon energy is coded in 
colour, ranging from 0.3 keV (red) to 5.0 keV (blue). As 
this image has not been corrected for exposure variations, 
the CCD boundaries are visible. The SNR N157B, the 
brightest source in the image, is saturated; an out-of-time 
event trail extends along the CCD column where the cen- 
tral neutron star is located. In the lower part of the image, 
sharp blue arcs appear. They are caused by single mirror 
reflections of photons from LMC X-l, which is outside the 
FOV. 



True images of the sky were obtained by projection of 
each photon onto the celestial plane, taking into account 
the different pointing directions and roll angles. The pro- 
jection takes care of the different spacing between the in- 
dividual CCDs as well as for the larger size of the pixels 
which arc furthest away from the readout nodes (Striider 
et al. |2001 ). The result of the projection, after removing 
events from bright pixels and a bright column, is shown in 
Fig. [t| The gaps between the CCDs show up in the image, 
as no exposure correction was applied. The exposure map 
(Fig.[|) shows that these areas also received some light. 
The exposure map was created in a similar way as the 
photon image, using the vignetting curve for 1.5 keV from 
the XMM Users' Manual (Dahlem [19991) . 



Apparent in Fig.[l] is a streak-like feature from the 
SNR N157B. This is an artefact of the operation of the 
CCD, which is sensitive to X-rays also during readout. In 
the full frame mode, 7% of all events are recorded dur- 
ing this period and get a displaced position along readout 
direction. This effect occurs everywhere on the detector, 
but out-of-time event trails are evident only for bright 
sources. Suppressing the trails does not only improve the 
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Fig. 2. Combined exposure map of all five observations. 
For the final image, only the inner part of the FOV was 
used, where the exposures overlap completely. Since the 
pointing position and roll angle of XMM-Newton were 
changed between the two sets of observations, almost all 
insensitive areas of the detector were exposed in the com- 
plementary observation. 

image cosmetically, but also reduces local variations of 
the background needed to apply the source detection al- 
gorithms. In the following we describe a straightforward 
method to suppress the contamination by out-of-time 
events. 

As the generation of out-of-time events is a random 
process, these events become uniformly distributed along 
each readout channel. The fraction of out-of-time events 
is given by the ratio of readout vs. frame time. Thus, if 
n is the total number of photons recorded in a particu- 
lar column, the contamination n CO nt of any pixel in this 
column by out-of-time events is 

^rcad 

^cont — T ' ™ 

Rrame ' ^pix 

where t re ad ; iframc and n p j X = 200 are the readout time, the 
frame time and the number of pixels per column. For the 
full frame mode, the values are t rea d = 5 ms and if ra mc = 
73 ms. Fig.|| shows the contamination of Fig. |l| by out- 
of-time events, calculated separately for the three colour 
bands in detector coordinates and then projected onto the 
sky, so that it can be subtracted from Fig.|l|. 

An additional source of undesired 'background' in 
Fig. |l| are sharp, blue arcs in the lower corner, which would 
systematically contaminate the spectra of the diffuse soft 
X-ray emission in this region. These are caused by sin- 
gle mirror reflections of photons from the bright X-ray 
binary LMC X-l, which is outside the FOV. In order to 




Fig. 3. Out-of-time event image of Fig.[y, constructed by 
distributing the events recorded in each readout column 
homogeneously across the column. The intensity of this 
image was artificially increased with respect to Fig. [j] by 
a factor of 20 for better visibility of these events. 
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Fig. 4. Single mirror reflections of LMC X-l in the lower- 
corner of Fig. |l|. Pixels affected during the first two obser- 
vations are marked in yellow, while pixels affected during 
the last three observations are coded in blue. The correc- 
tion was done individually for the two sets of observations. 



suppress these arcs, we identified the contaminated pixels 
from the hard images, where they are most distinct, and 
flagged them (Fig.|j). For the first light image, we replaced 
the flux of the affected pixels by interpolated values from 
the adjacent pixels. This approach seems to be justified, 
since the extent of one pixel (4'.' 13) is considerably smaller 
than the half energy width of the telescope (15'.'1), which 
implies that the fluxes of adjacent pixels cannot be inde- 
pendent of each other due to the point spread function of 
the telescope. 
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Fig. 5. Final version of the XMM-Newton EPIC pn first light X-ray colour image of the LMC, after applying all the 
data reductions described in the text. 



A final challenge for producing the first light image 
was the high dynamic range in intensity, which exceeds 
three orders of magnitude between the diffuse extended 
regions of low surface brightness and the point source in 
N157B. Since standard techniques like logarithmic inten- 
sity coding did not yield satisfying results, we designed a 
special code which preserves different shades of gray for 
the brighest sources. The final image, boxed with an equa- 
torial coordinate grid, is shown in Fig. ^[ 

3. Results and Discussion 

As Fig.|| shows, there is a wealth of different objects, 
point-like and extended up to 130 pc in size. The X-ray 



colours show large spectral differences from object to ob- 
ject. Based on this image, we accumulated photons from 
selected regions and performed spectral investigations of 
various sources, to utilize the full resolving power of the 
EPIC pn camera. We did this separately for the three 
data sets which were taken with the same observational 
setup (filter and pointing direction), because they require 
different detector response matrices. The best fit was ob- 
tained by modeling each of the three count rate spectra 
with the same candidate source spectrum but with the in- 
dividual normalization factor left as a free parameter (for 
details see Haberl et al. 



2001). Errors were determined at 



AX 2 



-2.3 with all parameters free. 
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3.1. 30 Dor C 

The largest object in the field is 30 Dor C, already known 
f rom radio observations and suggested by Mills et al. 
( |l984|) to be an SNR. In X-rays it appears as a com- 
plete ring of ~ 6 arcmin diameter. The size of the shell 
and its colour make this remnant the first of its kind out- 
side o ur Galaxy. In a separate paper (Aschenbach et al. 
2001a ) this unique object is investigated in detail. Close 
to the geometric center we detect a previously unknown 
point source, displayed in yellow, which emits mainly in 
the 0.3- 1.0 keV band. If this source is in the LMC, then 
its luminosity is of the order of 10 33 erg s _1 , and it might 
be the compact remnant of 30 Dor C, either a neutron 
star or a black hole. 



3.2. SN 1987 A 

Approximately 2 arcmin to the south-west of the shell 
border is the well-known SN 1987A. The XMM-Newton 
observations are very important with respect to both, the 
light curve and the X-ray spectrum, refining current mod- 
els of SN 1987A. Both aspects are presented and discussed 
in an accompanying paper by Aschenbach et al. ( |2001t; ). 



3.3. Active Galactic Nuclei behind the LMC 

A number of point sources have been discovered, displayed 
as blue dots in Fig.[|, the brightest of which is seen close 
to the image center. These sources have power law spectra 
and show strong low energy absorption. Both, the spectra 
and the fact that the central source is associated with a ra- 
dio point source, indicate that we see background sources, 
probably AGN, Nuclei, through and far beyond the LMC. 
A detailed analysis of these new X-ray source s is pr esented 
in an accompanying paper by Haberl et al. ( |200l| ) . 

3.4. N157B 



By far the brightest object in the field is the SNR N157B, 
located halfway from the image center to the edge in 
north-easterly direction. Afte r the discovery of th e 16 ms 
pulsar (Marshall et al. |1998| , Wang & Gotthelf |l998b| ) 
N157B has been classified as a Crab-like remnant with 
a central pulsar and a surrounding synchrotron nebula. 
In a preceeding paper Wang & Gotthelf ( 1998a ) showed 
that the spectrum is dominated by a power law with an 
indication of some thermal emission. The XMM-Newton 
spectrum of N157B (Fig. ^ demonstrates convincingly for 
the first time the presence of emission lines at low en- 
ergies, proving the presence of a thermal component. A 
model comprising two thermal components and a power 
law gives an acceptable fit to the data. Formally, ionization 
equilibri um m odels like Raymond-Smith (RS; Raymond 
& Smith 1977 ) or non equilibrium ionization models (NEI; 
Hamilton et al. |l983[) of XSPEC used for the thermal com- 
ponents at low energies cannot be discriminated. 




Energy (keV) 

Fig. 6. EPIC pn X-ray spectra of N157B for the three 
data sets with the same observational setup. The his- 
tograms represent the best fit model simultaneously ap- 
plied to the spectra. 

Fig. [6] shows the best fit two temperature RS mod- 
els with a power law model added. The temperatures are 
0.37t^5 keV and 0.12±0.01 keV, and the photon index 
of the power law is 2.83±0.03. With an assumed hydro- 
gen equivalent photoelectric absorption of 6-10 20 H-atoms 
cm~ 2 for the galactic foreground, we derive an additional 
H-equivalent column density of (1.94±0.04)-10 22 H-atoms 
cm -2 for the absorption in the LMC, taking elemental 
abundances of 0.5 solar into account. To fit the Mg emis- 
sion line at 1.3 keV clearly seen in the spectra, an over- 
abundance of 1.3 solar is required (for the cooler compo- 
nent). The NEI model yields similar values for the ab- 
sorption, the Mg abundance and the power law slope, but 
different temperatures (0.24 and 0.10 keV). 

The existence of the additional two-temperature ther- 
mal component suggests that N157B is not a Crab-like 
remnant but a composite remnant with an X-ray bright 
pulsar and synchrotron nebula but relatively faint ther- 
mal shell. The analysis of the thermal component, to be 
done in a separate paper, will be used to check for in- 
stance the consistency of the age derived from a Sedov- 
type approach with the pulsar characteristic spin-down 
age of P/(2P) ~ 5000 ys (Wang & Gotthelf |l998b|) . 

The power law component in the spectrum of N157B 
which has been discovered by Wang & Gotthelf ( 1998a ) 
using ROSAT PSPC and ASCA SIS data gives a pho- 
ton index of 2.3±0.1. The EPIC pn best fit value for the 
index is 2.83±0.03. N157B has been compared with the 
Crab. Both the Crab pulsar and the N157B 16 ms pul- 



sar (Marshall et al. 1998) show a flat power law spectrum 
with a photon index of 1.5 and 1.6±0.3, respectively. The 
spectrum of the Crab nebula is steeper by 0.5, generally 
believed to be due to synchrotron losses. If the spectrum 
of N157B would have a similar difference of ~ 0.5 between 
the pulsar and the nebula spectrum as suggested by the 
analysis of Wang & Gotthelf, the analogy with the Crab 
would be convincing. If on the other hand the EPIC pn 
power law slope is a better representation of the spectrum, 
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the analogy with the Crab is cleary violated and other ac- 
celeration and loss processes have to be considered. We 
note that the analysis of power law dominated spectra 
like in AGN show just minor discrepancies in the spectral 
index between ASCA and EPIC pn. Because of the steep 
spectrum of the nebula the physical conditions, either the 
electron injection spectrum, the electron transport speed 
or the magnetic field distribution, are likely to be different 
from those prevailing in the Crab. 

3.5. The Honeycomb nebula 

In Fig.^ X-ray emission from the Honeycomb nebula 
is seen 1 arcmin south-east of SN 1987A. The emission 
of the bright orange patch can be modeled by a ther- 
mal plasma with a temperature of 0.19±0.02 keV, su- 
perimposed on a power law of index 3.6±0.6. Increased 
Ne and Mg abundances (~2 times solar) are indicated. 
The hydrogen equivalent LMC absorption to this SNR 
is (2.2±0.9)-10 21 H-atoms cm -2 . Fig.0 shows the three 
EPIC pn spectra together with the best fit model. 




Energy (keV) 



Fig. 7. Same as Fig.O for the Honeycomb nebula. 



3.6. Interstellar medium in the LMC 

The north-western part of Fig. ^| appears remarkably dark 
and empty. But when we extract a spectrum from a 
44 arcmin 2 region north of 30 Dor C where no point 
sources and no out-of-time events from 30 Dor C affect 
the data, we clearly see thermal emission. The tempera- 
ture appears to differ from place to place over a range from 
~0.1 keV and ^0.3 keV, as was already concluded from 



ROSAT PSPC spectra (Sasaki et al. gOOlJ ). But for the 
first time we detect emission lines of oxygen, neon, mag- 
nesium and silicon (Fig.||). To fit the spectrum, we sub- 
tracted the background in a preliminary way, using source 
free areas from other EPIC pn observations. The result- 
ing spectrum is well represented by a 0.3l ; 2 keV RS 
type model with overabundances in O, Ne, Mg and Si of 
1.6±0.7, 5.4±1.5, 6.6±2.5 and 8.9±7.4 solar, respectively, 




Energy (keV) 

Fig. 8. EPIC pn X-ray spectrum of the faint, diffuse emis- 
sion in the north of Fig. ||, accumulated from the last two 
observations with thin filter. 

which is surprising in view of the generally low metallicity 
of the LMC. 
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